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ABSTRACT 

This  paper  reports  on  the  first  year’s  accomplishments  of  a  new  project  to  extend  explosion  models,  accounting  for 
source  medium  damage,  to  high  frequencies.  This  project  builds  on  a  previous  three -year  project  which  developed 
kinematic  source  descriptions  for  such  models  and  investigated  the  seismic  radiation  for  frequencies  below  ~0.2  Hz. 
Direct  and  indirect  effects  of  shock  waves  including  free  surface  interactions  cause  material  damage  which  in  general 
contributes  volumetric,  compensated  linear  vector  dipole  (CLVD),  and  double-couple  (DC)  sources  of  seismic 
radiation.  The  past  project  studied  two  important  features  of  these  models:  (i)  Rayleigh  waves  excited  by  a  CLVD 
source  and  the  impact  on  performance  of  m^-Ms  discrimination  and  (ii)  volumetric  moment  caused  by  damage  and  its 
impact  on  estimates  of  isotropic  moment  Mj  and  yield  estimation.  The  goal  of  the  current  project  is  to  advance  these 
models  with  a  sound  physical  basis  in  order  to  study  S  wave  generation  for  high  frequencies.  The  tasking  is  broken 
down  into  three  major  phases:  (1)  develop  and  quantify  mathematical  descriptions  of  “effective”  source  functions  for 
damage,  which,  up  until  now,  were  assumed  to  be  a  step  function,  (2)  test  and  validate  models  against  observed 
spectral  features  of  P/S  ratios  for  the  0.5-3  Hz  band,  and  (3)  test  and  validate  models  against  scaling  observations  of 
P/S  ratios  for  the  1-10  Hz  band.  This  year’s  work  involved  empirical  studies  quantifying  the  amount  of  volumetric 
moment  contributed  by  source  medium  damage.  Quantification  of  Mj  takes  the  form  Mj  -  M,  ■  Kx,  where  Mt  is  the 
classical  moment  due  to  cavity  formation,  Kis  a  measure  of  damage  tied  to  CLVD  moment  Mciv({,  and  exponent  x  is  > 
0  for  Pahute  Mesa  explosions  (Patton  and  Taylor,  201 1).  For  Mc\vd  -  0,  K  =  1,  and  M/  =  Mt.  For  a  CLVD  with 
permanent  extensional  deformation  along  the  vertical  axis,  ;V/r/lyy  >  0,  K  >  1,  and  Mj  >  Mt,  where  excess  volumetric 
moment  is  contributed  by  damage  processes  due  to  bulking  and  shear  dilation  of  the  source  medium.  The  largest 
explosions  on  Pahute  Mesa  have  K  values  of  1  or  slightly  less;  thus,  damage  does  not  radiated  much  energy  at  long- 
periods,  but  it  certainly  can  at  high  frequencies.  At  lower  yields,  K  gets  as  large  as  2.5  and  Mj  as  much  as  6  times  Mt. 
I  show  results  of  modeling  Ms  scaling  that  exploit  ground  truth  yields  in  order  to  place  constraints  on  the  value  of  x. 

Increased  understanding  of  the  explosion  source  was  applied  to  investigate  source  medium  properties  and  bias  of 
the  North  Korean  (NK)  test  site  using  yield  :  depth  of  burial  tradeoff  curves  and  { ni/,  -  Ms }  double  differences 
(Patton,  201 1).  The  results  favor  high  P  wave  speeds,  supporting  the  supposition  that  the  granite  medium  in  which 
NK  tests  were  conducted  is  strong  and  intact.  For  the  2009  NK  test,  revised  tradeoff  curves  accounting  for  better 
estimates  of  P  wave  speed  and  mu  bias  predict  yields  of  7-8  kt  for  a  burial  depth  of  550  m.  This  yield  and  burial  depth 
give  a  scaled  depth  of  -300  m/kt 1/3 .  An  emplacement  scenario  where  NK  tests  were  over-buried  in  competent  granite 
supports  previous  claims  that  damage  due  to  shock-wave  interactions  with  the  free  surface  was  suppressed,  leading  to 
poor  discrimination  performance  of  m^  :  Ms  (Patton  and  Taylor,  2008).  This  work  also  validated  our  understanding 
that  damage  does  not  radiate  energy  effectively  at  long -periods  for  large  Pahute  Mesa  explosions. 

Finally,  source  specifications  for  numerical  simulations  of  “effective”  S  wave  source  functions  are  being  developed. 
These  specifications  quantify  the  damage  distributed  in  depth  from  the  free  surface  to  the  depth  of  burial.  Effective 
source  functions  for  S  waves  will  be  computed  assuming  direct  generation  and  indirect  generation  through  Rg 
scattering.  A  review  of  the  specifications  will  be  presented  for  discussion  and  comment. 
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OBJECTIVE 

The  objective  is  to  develop  new  analytical  explosion  source  models  based  on  seismic  moment  tensor  theory  for 
further  improvement  and  advancement  of  regional  seismic  discrimination  and  yield  estimation  technologies.  Such 
technologies  rely  heavily  upon  the  source  information  contained  in  high-frequency  shear  (.S')  waves.  The  use  of  coda 
waves  following  regional  5  phases  to  estimate  explosion  yield  is  one  example  of  an  emerging  technology  offering 
great  promise  for  improved  nuclear  monitoring.  Unfortunately,  an  understanding  of  how  explosions  excite  5  waves  is 
quite  limited,  and  a  widening  gulf  between  theory  and  practice  undermines  our  confidence  to  monitor  broad  areas  at 
small  yields.  The  new  models  will  provide  a  physical  basis  for  explosion-generated  5  waves  and  theoretical  insights 
for  advancing  yield  estimation  and  discrimination  capabilities,  thereby  closing  the  gulf  between  theory  and  practice. 

RESEARCH  ACCOMPLISHED 


This  project  builds  upon  spherical  (monopole)  explosion  source  models  developed  in  the  1970’s.  An  important 
feature  of  those  models  is  the  theory  relating  seismic  amplitudes  of  P  waves  directly  to  yield,  depth  of  burial,  and 
material  properties  of  the  source  medium.  The  theory  draws  upon  empirical  yield  scaling  behaviors  of  key  model 
constructs,  such  as  the  elastic  radius,  and  the  analytical  nature  of  these  models  facilitated  their  use  since  they  were 
easy  to  implement,  and  as  such,  widely  applied  to  study  the  explosion  source.  Their  application  continues  to  this  day, 
but  with  the  recognition  that  a  spherical  point  source  is  inadequate  to  explain  5-wave  generation. 

The  model  under  development  consists  of  a  monopole  explosion  source  with  linear  superposition  of  tectonic  and 
material  damage  sources.  Material  damage  occurs  as  sudden  changes  in  the  source  medium’s  elastic  moduli.  To  date, 
our  research  has  focused  on  long-period  applications  where  all  source  components  are  safely  assumed  to  be 
coincident  in  time  and  space,  and  share  the  same  source-time  histories.  The  current  project  will  tackle  high-frequency 
applications  where  those  assumptions  are  no  longer  valid.  Before  doing  so,  certain  key  source  parameters  must  be 
quantified,  and  among  them  is  the  isotropic  or  volumetric  moment. 

There  are  two  sources  of  volumetric  moment  in  the  new  source  model,  one  from  cavity  formation  and  the  other  from 
damage  processes  (Patton  and  Taylor,  2011;  hereafter  PT1 1).  Since  tectonic  release  is  modeled  with  a  double-couple 
force  system,  it  does  not  contribute  volumetric  moment.  The  moment  tensor  for  a  monopole  explosion  source  is 


M  • 


1  0  0 
0  1  0 
0  0  1 


(1) 


where  M  is  a  volumetric  moment.  Assuming  an  incompressible  source  medium,  M  =  Mt,  the  moment  due  to  cavity 
formation,  where  Mt  =  per- Vc  (p  and  a  are  density  and  P  wave  speed  of  the  medium,  and  Vc  is  cavity  volume; 
Mueller  and  Murphy,  1971;  hereafter  MM71).  Vc  can  be  calculated  from  cavity  radius  scaling  models  (e.g..  Heard 
and  Ackerman,  1967;  Denny  and  Johnson,  1991;  hereafter  referred  to  as  HA67  and  DJ91).  The  moment  tensor  for 
damage  can  be  expressed  as  a  vertical  dipole  force  (Knopoff  and  Randall,  1970;  Ben-Zion  and  Ampuero,  2009). 


M, 


0  0  0 
0  0  0  . 
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where  Mtj  is  a  positive  number  (as  is  Mt  for  an  explosion),  z.M ,i  is  the  moment  of  the  dipole  force,  and  7  is  a  function 
of  a  new  source  parameter  K,  z  =  z(K).  Recall  that  K  is  defined  in  PT1 1  as 


2  M_ 


K  = 


MXX  +  Myy 


(3) 


where  Mxx ,  Mvv,  M 7Z  are  diagonal  elements  of  the  full  moment  tensor  source.  Decomposing  the  damage  source. 
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The  second  term  on  the  right  is  a  compensated  linear  vector  dipole  (CLVD)  with  moment  2zMJ2  and  vertical  axis  of 
symmetry  in  extension  if  z  >  0.  The  net  volumetric  moment  of  explosion  and  damage  sources  is  Mt  +  zMJ3.  The  ratio 
of  CLVD  moment  to  net  volumetric  moment  equals  2zM  J(~3M  t  +  zMj).  Assuming  tectonic  release  is  mainly  in  the 
horizontal  plane  and  contributes  little  to  ;V/ _ compared  to  explosion  and  damage  sources,  then 

M, 

z{K)  =  and  zMd  =  Mt  (K-V  -  (5) 

which  makes  use  of  equation  6  in  PT1 1 .  The  net  volumetric  moment  Mj  equals  Mt  •  (K+2)/3.  When  K  =  1,  the 
moment  of  the  dipole  source  is  zero,  and  Mj  =  M For  K>  1,  zM^  >  0,  and  M/  >  Mt\  meanwhile  if  K  <  1,  z.Mtj  <  0,  and 
Mj  <  Mt.  Thus,  M 1  exceeds  the  moment  due  to  cavity  formation  if  the  value  of  source  parameter  K  is  greater  than  one. 

Predictions  of  Mt  were  made  by  PT1 1  for  comparison  with  Mj  estimates  of  Pahute  Mesa  explosions.  The  results 
showed  that  Mj  >  Mt  for  explosions  with  K  >  1,  consistent  with  an  additional  source  of  volumetric  moment  due  to 
damage.  To  quantify  how  much  additional  moment  was  present,  PT1 1  adopted  a  power-law  model  of  the  form,  Mj  = 
M,  ■  Kx,  and  attempted  to  estimate  the  exponent  x.  Unfortunately,  Mj/Mt  ratio  observations  show  too  much  scatter  to 
allow  a  very  reliable  determination  of  x.  Due  to  the  importance  of  quantifying  the  excess  moment  and  testing 
observations  against  the  model,  Mj=  M ,  •  (A'+2)/3,  additional  work  has  been  carried  out  using  Ms  yield  scaling 
observations  for  Pahute  Mesa  explosions  with  ground-truth  yields  to  constrain  the  range  of  acceptable  x  values  more 
tightly  than  Mj/Mt  observations  were  able  to. 

Constraints  on  x  from  Ms  Yield  Scaling  Observations 

The  Rayleigh-wave  reduced  excitation  spectrum  A(co)  is  a  measure  of  the  azimuth-independent  radiation  of  the 
seismic  source  (Patton,  1988).  It  is  extracted  from  a  linear  decomposition  of  path-corrected  complex  amplitudes  (real 
and  imaginary  parts)  for  a  station  network  into  five  constituent  azimuthal  radiation  patterns  (1,  cos0,  sin0,  cos20, 
sin20;  Romanowicz,  1982).  The  A  spectrum  is  related  to  the  strength  and  polarity  of  the  constant  term.  For  long- 
periods  T  (T  =  27t/co,  to  is  angular  frequency),  the  A  spectrum  can  be  written  in  terms  of  Rayleigh-wave  excitation  and 
source  parameters  for  the  model  under  development  (Patton  and  Taylor  2008;  hereafter  PT08) 


A(a>)  =  f  (K)  ■  M j  ■  G^co)  , 


(6) 


where  G  |  ( &))  is  a  Green’s  function  and /OK)  =  (6-2K)/(2+K).  Using  the  DJ91  scaling  model,  Mt  can  be  expressed  as 
..  .  0.21  0.85  ,-0.79 

Mt  =  A  •  p  •  a  •  /z  •  W  ,  (7) 


where  h  is  depth  of  burial,  and  A  is  a  constant  of  proportionality.  Using  the  power  law  model  Mj  =  Mt  ■  Kx ,  then 
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A(co)  =  A  ■  f(K)  ■  Kx  ■  Cl  ■  W  ■  Gj(CO)  ,  (8) 

where  Cl  =  Cl(h)  =  p°'~ 1  •  a°'85  •  h  °'7  >  and  p  =  p(/z),  a  =  a (h).  PT08  developed  an  empirical  template  A/co)  for  A 
spectra  using  an  ensemble  of  Pahute  Mesa  explosions.  The  gain  A'(.  is  a  multiplicative  scaling  factor  applied  to  the 
template  to  match  the  A-spectrum  amplitude  of  the  zth  explosion.  Between  0.1  and  -0.025  Hz,  the  template  matches 
the  frequency  dependence  of  G  |  computed  for  a  Pahute  Mesa  velocity  model.  The  model  predicts  little  variation  of 
G'|  for  burial  depths  between  500  and  1500  m.  Thus,  Af(  to)  ~  y  ■  G  |  ( 0J)  and  A,.(  to)  ®  A'  -  ■  y  ■  G  ( ( 0J ) ,  where  y  is  an 
event-independent  scaling  factor  with  physical  units  of  moment  in  Nm.  Substituting  for  A(co)  into  equation  (8),  the 
gain  factor  A'  •  for  the  zth  explosion  is 

A'f»  A' ■/(*,.)  ■K^CliWi  ,  (9) 

where  A'  (=  A/y)  is  independent  of  yield.  I  will  establish  empirically  that  log[  A'  ]  and  network  Ms  scale  one-to-one 
for  nuclear  explosions  detonated  on  Pahute  Mesa.  This  scaling  equivalence  will  allow  A7v-yield  scaling  to  be 
interpreted  with  the  new  explosion  source  model  since 

Msi~C,log[Wi]~log[f(Ki)]  +  xlog[Ki]  +  log[Q.i]  +  log[W(]  ,  (10) 

where  C,  is  the  apparent  yield  scaling  exponent  observed  for  Ms.  Density  and  P  wave  speed  of  the  source  medium  (and 
hence  Cl)  vary  with  burial  depth,  and  therefore  with  yield  for  a  containment  practice,  such  as  hl  ~  120 -W)  .  Since  K 

shows  a  systematic  dependence  on  yield  (PT1 1),  all  terms  on  the  right-hand  side  of  equation  (10)  will  be  expressed 
with  a  power-law  dependence  as  follows 


K{  ~  W“  , 

ai) 

6  -  2  K{  b 

/(*■«)  -  27 ir.~wi  -and 

(12) 

^  0.21  0.85  ,-0.79  c 

ai  =  Pi  ■  ai  ■  hi  ~  Wi  ’ 

(13) 

where  exponents  a,  b,  c  and  C  will  be  estimated  from  measurements  or,  in  the  case  of  Qr  from 
cube-root  containment  rule.  The  yield  exponent  C,  is  related  to  the  other  exponents  as  follows 

a  velocity  model  and 

C,  =  ax  +  b  +  c  +  1  , 

(14) 

and  solving  for  x 

v  _  (C-l )-{b  +  c) 
a 

(15) 

Values  of  x  consistent  with  observed  yield  scaling  of  Ms  for  Pahute  Mesa  explosions  will  be  found,  but  first  the 
scaling  equivalence  of  network  Ms  and  log  [A'  ]  must  be  shown. 

Scaling  Equivalence  of  Ms  and  log[A'  ]  and  Yield  Scaling  Analysis 

Figure  1  displays  plots  of  log[  A’  ]  against  Ms  for  five  different  published  data  sets.  All  network  Ms  determinations, 
with  the  exception  of  Ms  reported  by  Selby  et  al.  (201 1),  scale  one-to-one  with  log  [A1  ].  An  average  network  Ms  was 
computed  using  Jih  and  Baumstark  (1994)  Ms  values  to  normalize  the  baseline  offsets  seen  in  Figure  1  (a-d).  The 
average  Ms  plotted  against  log[A'  ]  has  a  slope  of  1.00  and  a  l-o  formal  error  from  regression  analysis  of  ±  0.03. 
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56  Pahute  Mesa  Explosions  46  Pahute  Meso  Explosions 


Figure  1.  Plots  of  various  network  Ms  data  sets  against  log  [A']  for  Pahute  Mesa  explosions.  Ms  data 
sets  are  taken  from  (a)  Stevens  and  Murphy  (2001),  (b)  Woods  and  Harkrider  (1995),  (c) 

Jih  and  Baumstark  (1994),  (d)  Bonner  et  al.  (2003),  and  (f)  Selby  et  al.  (2011).  Average 
network  Ms  are  based  on  the  first  four  data  sets  with  a  baseline  established  arbitrarily  by 
Jih  and  Baumstark.  Open  circles  in  (e)  are  explosions  with  only  one  data  set  reporting  a 
network  Ms;  closed  circles  have  multiple  data  sets  reporting. 

Note  a  change  in  vertical  scale  in  Figure  le,  and  the  arbitrary  shift  in  the  line  with  a  slope  of  1  to  facilitate  comparison 
to  the  data.  Average  Ms  are  available  for  76  Pahute  Mesa  explosions.  In  addition,  Ms  values  were  available  for  56 
explosions  by  converting  log[  A'  ]  measurements  to  equivalent  Ms  values  by  adding  a  constant  offset  of  4.655 
magnitude  units  (mu)  determined  from  the  data  plotted  in  Figure  le. 
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In  order  to  mitigate  effects  of  the  water  table  on  the  results,  the  yield-scaling  analysis  was  restricted  to  explosions 
below  the  nominal  standing  water  level  on  Pahute  Mesa  (640  m).  Thirty-four  explosions  meet  this  restriction,  Kash 
with  the  shallowest  burial  depth  (645  m)  and  Muenster  the  deepest  (1452  m).  Rex,  a  high-leverage,  overburied  shot 
with  an  announced  yield  of  19  kt  and  burial  depth  of  671  (Springer  et  al.,  2002;  -250  m/kt  ),  was  omitted  from  the 
analysis,  bringing  the  total  down  to  33.  Of  those  33  shots,  27  have  K  values  and  25  have  inferred  Ms  from  log[  A'  ]. 


Three  coupling  scenarios  are  considered:  (1)  no  coupling  variations  with 
burial  depth  (DOB;  uniform  coupling),  (2)  coupling  appropriate  for  P 
waves  based  on  assessments  of  Murphy  (1996),  and  (3)  coupling  based 
on  a  sigmoid  curve  fit  E(/z)  to  surface  wave  observations  (PT1 1). 

According  to  Murphy  (1996,  p.  225),  “...for  explosions  in  dry,  porous 
media,  such  as  dry  alluvium  or  tuff,  ...  the  average  mb  values  for  a  given 
yield  are  lower  than  those  in  hard  rock  by  about  0.50  ±  0.25  mu.” 

Coupling  variations  for  P  waves  were  based  on  the  same  sigmoid  curve 
used  for  surface  waves,  but  adopting  0.5  mu  for  the  shallow-depth 
asymptote  (coupling  factor  of  10_0  5  =  0.32).  The  shallow-depth 
asymptote  of  £(/z)  has  a  coupling  factor  of  -0.04.  Coupling  curves  for  P 
and  surface  waves,  plotted  in  Figure  2,  differ  by  about  a  factor  of  two  at  a  burial  depth  of  640  m.  P  wave  coupling 
factors  weight  the  down-going  rays  more  heavily  than  do  surface  waves  which  require  a  free  surface  for  their 
existence.  Down-going  rays  encounters  water-saturated  media  where  coupling  improves,  while  surface  waves  sample 
the  integrated  effects  of  rays  leaving  at  all  take-off  angles.  As  such,  one  might  expect  surface  wave  coupling  factors  to 
be  smaller  than  their  P  wave  counterparts  for  Nevada  National  Security  Site  (NNSS,  formerly  the  Nevada  Test  Site 
NTS)  media. 


Figure  2.  Coupling  factors  for  P  waves 

(black)  and  Rayleigh  waves  (red). 


Estimates  of  model  parameters  for  linear  regressions  involving  official  yields  are  provided  in  Tables  1  through  3. 
Table  1  lists  estimates  of  the  yield  scaling  exponent  C,.  Results  are  provided  for  the  combined  data  sets  of  average  Ms 
and  inferred  Ms,  a  total  of  58  measurements,  as  well  as  data  sets  of  just  33  average  Ms  measurements  and  25  inferred 
Ms  from  log[  A'  ].  While  the  results  are  consistent  across  all  three  data  sets,  systematic  variations  in  estimates  of  C  for 
the  different  coupling  scenarios  are  found  as  expected.  The  scaling  exponents  for  surface  wave  coupling  are  smaller 
than  their  counterparts  for  P  wave  coupling  which  in  turn  are  smaller  than  for  uniform  coupling. 


Table  1:  Mv-log[  W\  Scaling  for  Three  Coupling  Scenarios,  Pahute  Mesa 


Data  set 

Uniform 

P  wave 

Surface  wave 

58  measurements 

33  with  Ms 

25  with  log[A'] 

£  =  1.144  ±0.036 
o  =  0.111 

CC  =  0.97 

C  =  1.039  ±0.035 

0  =  0.107 

CC  =  0.97 

C,  =  0.853  ±  0.038 

0  =  0.118 

CC  =  0.95 

33  shots  with  Ms 
h  >  640  m 

£  =  1.143  ±0.052 

0  =  0.123 

£  =  1.039  ±0.049 

0  =  0.115 

t,  =  0.850  ±  0.051 

0  =  0.120 

25  shots  with  log[A'] 
h  >  640  m 

£  =  1.148  ±0.049 
o  =  0.098 

£=  1.046  ±0.050 
o  =  0.100 

C,  =  0.860  ±  0.059 

0  =  0.118 

Table  2  summarizes  scaling  results  for  K,f(K ),  and  Ll(li).  The  dependence  of  iHh)  on  yield  was  obtained  for  a  cube- 
root  containment  practice  using  layered  models  based  on  studies  of  Earth  structure  on  Pahute  Mesa  by  Ferguson  et  al 
(1994)  and  Leonard  and  Johnson  (1987).  Estimates  of  density  and  P-wave  speed  appropriate  for  surface  wave 
excitation  were  obtained  using  a  method  that  averages  material  properties  over  the  elastic  radius  (see  PT1 1). 
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Table  2:  Yield  Regressions  on  K,f{K),  and  il(h) 


K 

m 

Q.(h) 

n  =  27,  h  >  640  m 
a  =  -0.291  ±0.044 

a=  0.092 

CC  =  0.80 

n  =  27,  h  >  640  m 
b  =  0.262  ±  0.046 
cr=  0.096 

CC  =  0.76 

n  =  33,  h  >  640  m 
c  =  -0.151  ±0.011 

a  =  0.025 

CC  =  0.93 

n  =  number  of  explosions;  a,b,c  =  scaling  slopes;  a  =  standard  deviation 
of  residuals;  CC  =  data  correlation 

With  estimates  of  scaling  slopes  a ,  b,  and  c  so  obtained,  equation  (15)  can  be  used  to  compute  values  of  x,  the 

exponent  on  K.  The  results  are  provided  in  Table  3.  While  x  does  not  show  variation  between  different  data  sets  due  to 
consistent  M,-yield  scaling  rates  £  in  Table  1,  there  are  significant  variations  for  the  three  coupling  scenarios. 
Uniform  coupling,  which  we  know  is  not  appropriate  for  any  test  area  at  NNSS,  has  a  small  and  slightly  negative 
exponent.  Coupling  scenarios  based  on  P  wave  and  surface  wave  estimates  (Figure  2)  have  positive  exponents 
consistent  with  surplus  volumetric  moment  above  and  beyond  what  cavity  formation  predicts.  Based  on  surface  wave 
coupling,  Ms  yield  scaling  favors  an  exponent  of  ~1.  If  surface  wave  coupling  factors  are  biased  for  some  reason,  the 
exponent  might  be  as  small  as  0.2  based  on  P  wave  coupling,  but  such  coupling  for  surface  waves  seems  unlikely  for 
reasons  stated  above.  Thus,  Ms  scaling  results  presented  herein  bound  the  exponent  on  K  to  be  greater  than  0.2  but 
perhaps  less  than  0.9.  This  range  is  significantly  smaller  than  that  suggested  by  Mj/Mt  observations  (x  >  1 ;  PT1 1). 


Table  3:  Estimation  of  Scaling  Exponent  x 


Uniform 

P  wave 

Surface  wave 

Data  set 

c 

X 

c 

X 

c 

X 

58  measurements 

1.14 

-0.1 

1.04 

+0.2 

0.85 

+0.9 

33  shots  with  M s 

1.14 

-0.1 

1.04 

+0.2 

0.85 

+0.9 

25  shots  with  log  [A'] 

1.15 

-0.1 

1.05 

+0.2 

0.86 

+0.9 

x  =  [  (£  -  1)  -  (b  +  c)  ]  /  a;  estimates  of  a ,  b. 

and  c  provided  in  Table  3 

As  shown  in  the  previous  section,  the  source  model  predicts  the  following  relationship  between  net  volumetric 
moment  Mj  and  moment  due  to  cavity  formation  Mt 

For  1  <  K  <  2.5,  where  2.5  is  the  upper  limit  measured  for  Pahute  Mesa  explosions,  the  function  (K  +  2)/3  is  best  fit  by 
the  power-law  model,  0 .977-/f°'453.  Thus  bounds  on  x  (0.2  <  x  <  0.9)  estimated  from  Ms  yield  scaling  observations  are 
consistent  with  the  source  model  prediction  (0.453).  K  values  of  1.5,  2.0,  and  2.5  predict  net  volumetric  moment  Mj 
exceeds  cavity  formation  moment  Mt  by  -17,  33,  and  50%,  respectively,  due  to  damage  of  the  source  medium.  A 
back-of-the -envelope  calculation  for  a  100-kt  explosion  with  K  value  of  2  translates  a  33%  moment  increase  into 
-0.1%  volumetric  dilation  of  damaged  source  medium  confined  to  an  inverted  cone  extending  from  the  cavity  formed 
by  the  explosion  to  the  free  surface  and  an  circular  area  around  ground  zero  where  spallation  is  strongest  (-160  scaled 
meters). 
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Inferences  of  Material  Properties  and  mb  Bias  of  the  North  Korean  Test  Site. 

Even  with  complete  understanding  of  explosion  source  mechanisms,  accurate  yield  estimation  still  requires  reliable 
knowledge  of  source  medium  velocity  structure  and  near-source  propagation  effects.  For  new  test  sites  in  remote, 
seismic-inactive  areas,  such  knowledge  can  be  difficult  to  obtain,  requiring  extensive  data  collection  efforts  and 
painstaking  observational  work  that  can  take  years  before  producing  reliable  results.  The  North  Korean  (NK)  test  site 
is  a  prime  example  of  just  this  scenario  where  little  is  known  about  the  material  properties  of  the  source  region  and 
the  efficiency  of  wave  propagation  under  the  northern  Korean  Peninsula. 

In  this  section,  I  will  show  how  source  theory  can  be  exploited  to  make  inferences  about  source  medium  wave  speeds 
and  nib  bias  for  the  NK  test  site  using  (1)  yield  :  depth-of-burial  tradeoff  curves  (W:h  TOCs)  and  (2)  {mb-Ms)  double 
differences  between  the  2009  NK  test  and  explosions  conducted  at  NNSS.  Herein  W:h  TOCs  based  on  measurements 
of  mb  and  Mj  are  presented  for  the  2009  test.  A  fundamental  assumption  made  to  derive  TOCs  is  that  the  test  can  be 
approximated  by  a  pure  explosion  source.  Measurements  of  mb,  Ms ,  and  the  results  of  moment  tensor  inversions 
suggest  that  the  sources  of  both  NK  tests  are  nearly  pure  explosions  with  relatively  minor  energy  radiated  at  long 
periods  from  tectonic  release  and  source  medium  damage  (e.g..  Ford  et  al.,  2009;  Shin  et  ah,  2010;  PT08,  PT11). 
Sharpe’s  (1942)  theory  will  be  used  to  relate  seismic  radiation  to  the  reduced  displacement  potential  t|/(x)  where  T  is 
reduced  time  t  -  r/a,  r  is  radial  distance,  and  a  is  the  speed  of  P  waves  in  the  medium.  I  make  use  of  two  well-known 
scaling  models,  one  due  to  HA67  (used  by  MM71)  and  the  other  due  to  DJ91.  These  models  predict  different  TOCs. 
The  material  presented  in  this  section  is  abstracted  from  Patton  (2011). 

Figure  3  shows  W:h  TOCs  for  the  2009  NK  test  based  on  measurements 
of  mb  and  Mb  from  Selby  et  ah  (201 1).  Assumed  material  property  values 
for  the  NK  test  site  are  those  adopted  by  Koper  et  ah  (2008;  hereafter 
KHB08):  pNK  =  2500  kg/m3,  aNK  =  5100  m/s,  vNK  (Poisson  ratio)  = 

0.2354,  and  GP^k  (gas  porosity)  =  0.5%.  The  test  site  bias  SXK  was 
assumed  to  be  zero.  Of  interest  are  the  offsets  between  TOC (mb)  and 
TOC  {Mj)  for  both  scaling  models,  where  TOC  (Mb)  >  TOC  (mb).  Not  only 
are  the  offsets  significant,  but  so  are  the  predictions  of  the  depth 
dependence  of  yield.  Important  questions  to  answer  are:  (1)  what  causes 
the  offsets,  and  (2)  which  set  of  TOCs,  HA67  or  DJ91,  is  a  better 
representation  of  the  truth.  This  study  addresses  the  first  question.  I 
assert  that  the  offsets  are  caused  by  systematic  errors  due  to  inaccurate 
medium  parameters  and  mb  bias  adopted  for  the  computation  of  TOCs  in 
Figure  3.  Measurement  errors  in  mb  and  Mj  are  assumed  to  be  second 
order  compared  to  systematic  errors.  As  such,  offsets  provide  constraints 
on  medium  properties  and  SNK.  Two  constraint  equations  involving  aXK, 
vNK,  and  SNK  will  be  derived  from  the  TOCs  and  from  {mb- Ms }  double 
differences  between  the  2009  NK  test  and  calibration  explosions 
detonated  on  Pahute  Mesa.  {mb  -  Ms}  double  differences  between  the  2009  test  and  NNSS  explosions  constrain  a 
three-dimensional  surface  relating  free  parameters  aNK,  vNK,  and  8XK.  A  fourth  dimension  of  parameter  space,  pNK, 
was  eliminated  by  the  use  of  a  p-a  relationship.  NNSS  serves  as  a  reference  test  site  where  material  properties 
aNNSS’  VNNSS  ani3  test  site  bias  SNNss  are  calibrated,  and  the  explosion  sources  are  well-quantified.  NNSS  sources 
must  also  be  pure  explosions  to  the  first  degree,  as  the  NK  explosions  are.  Two  selection  criteria  were  used:  (1)  mb  > 
5.8  and  (2)  h  >  -750  m.  PT1 1  found  that  large  Pahute  Mesa  explosions  meeting  these  criteria  have  measured  Mj  that 
are  consistent  with  the  predictions  of  classical  moment  Mt  based  on  cavity  formation  alone.  Twenty  explosions  met 
these  criteria. 


Figure  3.  W:h  TOCs  for  the  2009  NK  test. 

Two  sets  are  shown:  black  curves 
based  on  HA67  scaling  and  red 
curves  for  DJ91  scaling.  Solid 
lines  are  TOC(M/)s,  dotted  lines 
are  TOC(/nj).  *  symbols  indicate 
TOC  obtained  by  Murphy  et  al. 
(2010)  from  modeling  regional  P 
wave  source  spectra. 
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An  expression  for  { nij,  -  Ms]  double  differences  between  tests  conducted  at  the  NK  test  site  and  NNSS  is 
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(see  Patton,  201 1).  \)/NK  1  is  shorthand  for  ancl  (  =  (p/a)2.  tOj  and  0)o  are  frequencies  where  ni/,  and  Ms  are 

measured.  Equation  (17)  involves  m f,  biases  for  both  test  sites.  S^nss  was  estimated  from  nij/P)  -  nij/Lg)  residuals 
for  the  20  calibration  shots  (-0.27  ±  0.02  mu).  Ratios  \|/NNSS  i  /  Vnnss  1  were  measured  on  potentials  computed  with 
the  MM71  model.  Since  HA67  and  DJ91  predict  similar  cavity  radii  for  standard  burial  practice,  the  ratios  were  not 
affected  by  the  choice  of  scaling  model.  Ratios  \|fNK  |  /  \|/NK  2  were  measured  for  six  realizations  of  MM7 1  spectra  for 
granite  explosions  with  yields  3,  6,  and  12  kt  and  two  NK  velocity  models,  one  with  aNK  of  5100  m/s  and  one  with 
5600  m/s.  The  average  log  ratios  for  the  two  velocity  models  are  0.10  and  0.06  mu,  respectively.  For  NK  models, 
spectra  were  computed  for  a  scaled  burial  depth  of  250  m/kt  .  As  such,  the  DJ91  model  predicts  smaller  cavity  radii 
than  HA67,  hence  higher  corner  frequencies  and  smaller  ratios.  I  adopted  a  value  of  0.06  mu  for  log[\|/NK  l  /  V|/NK  2], 
which  is  0.1  mu  smaller  than  log[\(/NNss  l  /  ¥nnss  li-  imb  -/WV}NK  equals  1.04  mu  and  the  mean  {m^  -MS}NN 35  for 
20  Pahute  Mesa  explosions  is  1.19  ±  0.04  mu.  After  making  substitutions,  equation  (17)  reduces  to 
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Using  appropriate  values  for  Cnnss  anc*  aNNSS’  Cnnss  '  aNNSS  =  7.37xl06  .  Solving  for  aNK,  a  three-dimensional 
surface  is  obtained 
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Figure  4.  Contours  (black  dotted  lines) 

satisfying  the  constraint  TOC(M/)  = 
TOC(/nfi).  Contours  from  equation 
(19)  are  plotted  for  selected  8NK 
ranging  from  -0.05  to  0.25  mu  in 
steps  of  0.05  (thin  black,  heavy 
black,  blue,...,  red).  A  plausible 
upper  bound  on  aNK  for  near¬ 
surface  granite  is  shown  by  the 
horizontal  dash  line  (6000  m/s). 


An  offset  Ow  is  defined  as  a  multiplicative  factor  such  that  TOC (Mf)  = 
Ow  ■  TOC (mb).  Equations  for  offsets,  (/[(/(HA 67)  and  Owi  D.I9 1 ),  lead 
to  expressions  for  tX\jK  which  exhibit  simple  inverse  relationships 
where  aNK  increases  as  decreases  or  vice-versa. 
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a^K( DJ9 1 )  is  5.5%  faster  than  a^K(HA67).  Using  equation  (19),  the 
formulas  for  aNK  can  be  written  in  terms  of  (NK 
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(21) 


These  formulas  are  plotted  in  Figure  4  along  with  contours  of  constant 
;«/,  bias  from  equation  (19).  P  wave  speeds  at  the  NK  test  site  inferred 
from  both  scaling  models  are  considerably  faster  than  5100  m/s  assumed  by  KHB08  if  8\jK  is  small  (<  0.15  mu),  or 
under  the  reasonable  assumption  that  vNK  for  intact  granite  is  less  than  0.3. 
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Confirmation  that  revised  TOCs  are  completely  consistent  for  the 
solutions  in  equations  (20)  and  (21)  is  provided  in  Figure  5.  Now 
the  TOCs  overlay.  These  TOCs  bring  out  several  important  points 
about  the  solutions.  The  first  underscores  what  the  contours  in 
Figure  4  revealed,  that  source  medium  P  wave  speeds  must  be  high 
when  SNK  is  small.  Since  compliance  of  the  medium  stiffens  as 
aNK  and  pNK  increase,  the  elastic  response  is  reduced  giving  a 
smaller  P  wave  amplitude  for  a  given  applied  pressure.  Thus  the 
TOCs  for  zero  bias  predict  larger  yields  than  initial  TOC(m^)s  did 
(compare  dotted  curves  in  Figure  3  with  the  curves  for  8NK  =  0  in 
Figure  5).  The  second  point  is  that  increasing  8NK  lowers  the  yield 
predicted  by  TOCs.  A  8NK  of  0.3  mu  would  increase  the  yield  by 
about  a  factor  of  two  for  DJ91  scaling  were  it  not  for  compensating 
effects  of  a  more  compliant  source  medium  and  larger  cavity  radius 
since  aNK  decreases  as  SNK  increases.  For  DJ91  scaling,  TOCs 
decrease  a  factor  of  1.4  for  a  0.3  mu  increase  in  8NK  (Figure  5). 


Figure  5.  W:h  TOCs  for  a  range  of  8NK  values. 

TOC(M/)  and  TOC(wft),  plotted  with 
solid  and  dotted  lines,  overlay  as  they 
should.  5nk  values  are:  black,  blue, 
cyan,  magenta,  and  red  for  0.0,  0.05 
(HA67  only),  0.1,  0.2,  and  0.3  mu 
(DJ91  only),  respectively. 


Figure  6.  Revised  TOCs  for  2009  NK  test. 
TOC  for  HA67  scaling  (black) 
compares  well  with  Murphy  et 
al.’s  TOC  shown  by  *  symbols  on 
the  graph,  and  final  yield,  DOB 
shown  by  symbol.  Plotted  in  red 
are  “end  member”  solutions  for 
DJ91  scaling  with  minimum  and 
maximum  SNK  (solid  and  dash 
lines).  Magenta  lines  are  for 
constant  scaled  depth  of  burial 
(e.g.  120, 169.7,  240,  339.4  m/kt1/3). 


Revised  W:li  TOCs  are  plotted  in  Figure  6  for  the  solutions  indicated 
with  large  black  dots  in  Figure  4.  TOCs  computed  using  FIA67  scaling 
for  8nk,  0CnK’  vnk  °f  0.034  mu,  5500  m/s,  and  0.174  are  consistent 
with  the  TOC  of  Murphy  et  al.  (2010).  Final  estimates  of  yield  and 
depth  of  burial  reported  by  Murphy  et  al.  are  4.6  kt  and  550  m.  Two 
solutions  for  DJ91  scaling,  one  for  zero  bias  and  one  for  what  I  believe 
is  an  upper  bound,  0.15  mu,  are  plotted.  These  solutions  have  aNK  of 
5910  and  5450  m/s,  and  vNK  of  0.178  and  0.323,  respectively,  and 
bracket  a  plausible  range  of  TOCs.  Lines  of  constant  scaled  depth  of 
burial  are  also  plotted,  showing  graphically  that  scaled  burial  depths 
vary  along  the  TOCs  from  -80  m/kt1/3  to  more  than  500  m/kt1/3  for  the 
range  of  burial  depths  plotted.  It  is  apparent  from  these  results  that  the 
choice  of  scaling  model  is  quite  important  for  accurate  yield  estimation 
once  explosions  become  over-buried.  Hydrodynamic  simulations  for 
competent  granite  from  the  study  of  Rougier  et  al.  (201 1)  favor  DJ91 
scaling  and  provide  estimates  of  a  lower  bound  on  yield  and  depth  of 
burial.  For  the  2009  NK  test,  these  estimates  are  5  to  6  kt  and  -380  m. 
At  550  m,  the  depth  favored  by  Murphy  et  al.  (2010),  the  yield  from 
this  study  is  7  to  8  kt. 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  project  advances  new  explosion  source  models  by  accounting  for  the  effects  of  source  medium  damage  on 
radiated  P,  S  and  Rayleigh  wave  fields.  Volumetric  moment  due  to  damage  has  been  quantified.  Further  research  is 
needed  to  test  models  against  observed  P/S  ratio  spectral  modulations  in  the  0.5-3  Hz  band  and  high-frequency  yield 
scaling.  To  do  so,  the  next  step  will  build  suitable  time  histories  and  depth  distributions  for  damage  moment-rate 
functions.  Both  time  history  and  moment  will  depend  on  depth.  A  method  exploiting  new  insights  into  the  explosion 
source  led  to  constraints  on  material  properties  for  the  NK  test  site  and  estimates  of  yield  for  the  2009  NK  test. 
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